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Spectral Properties of Magnetic Chain Graphs 


P. Exner and S. Manko 


Abstract. We discuss spectral properties of a charged quantum particle 
confined to a chain graph consisting of an infinite array of rings under 
influence of a magnetic field assuming a (5-coupling at the points where 
the rings touch. We start with the situation when the system has a 
translational symmetry and analyze spectral consequences of perturba¬ 
tions of various kind, such as a local change of the magnetic field, of the 
coupling constant, or of a ring circumference. A particular attention is 
paid to weak perturbations, both local and periodic; for the latter we 
prove a version of Saxon-Hutner conjecture. 
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1. Introduction 

Quantum graphs are a subject both mathematically rich and useful in ap¬ 
plications; we refer to the recent monograph |BK13) for a broad panorama 
of this theory. One class of interest in this category are chain graphs hav¬ 
ing the form of an array of rings to which a quantum particle is confined. 
At the touching points of rings the wave functions are connected by suit¬ 
able matching conditions, one of the simplest possibilities being the so-called 
d-conditions described below. 

As long as such a chain is periodic, the spectrum of the corresponding 
Hamiltonian has a band structure, however, some of the band may be flat 
since the unique continuation principle of the usual Schrodinger theory in 
general does not hold for operators on graphs with a nontrivial topology. 
The picture changes when the periodicity is lost, local perturbation may give 
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rise to eigenvalues in the gap of the unperturbed spectrum associated with 
localized states, see e.g. |DET08] . 

The picture becomes even more intriguing if the particle in question 
is charged and exposed to a magnetic field, even in the simplest situation 
when the chain graph E is planar, the field is perpendicular to the graph 
plane, and no other potential is present. The Hamiltonian then acts as a one¬ 
dimensional magnetic Schrbdinger operator on each edge. Its domain consists 
of all functions from the Sobolev space and we assume that they obey 

the i5-coupling at the graph vertices which are characterized by conditions 

n 

■0^(0) = V'i(O) =:'0(0), ij = ^ ^'0»(O) = a'0(0), (1.1) 

i=l 

in a vertex where n edges meet. In our case we will have n = 4; the meaning 
of the quasiderivative ^ in the formula will be explained in the next section. 

In our previous paper |EM15) we discussed properties of a straight mag¬ 
netic chain and eigenvalues coming from local changes of the coupling con¬ 
stant a in (El). Our aim here is to analyze a broader class of perturbations 
including local variations of the magnetic field and edge lengths. We shall 
also discuss weak perturbation, not necessarily of local character. If those are 
periodic we are going to prove an analogue of the well-known Saxon-Hutner 
conjecture [AOHH051 Sec. III.2.3]. 

Let us briefly describe the contents of the paper. As we have mentioned, 
the periodic magnetic chain was analyzed in [EM 1 5] and we recall here the 
results only to the degree needed to make the present paper self-contained. 
Our main technical tool will be a conversion of the problem into a difference 
equation. Such a duality trick is well known |Ca97[ lEx971 IPaI3] . however, we 
have to work it out for our purposes; this will be done in Sec. El The results 
are contained in the next two sections. Eirst we show how the discrete spec¬ 
trum coming from compactly supported perturbations of various sorts can 
be found and discuss in detail several examples. Next, in Sec. O we analyze 
the weak-perturbation situation, first for compactly supported variations of 
the magnetic field and coupling constants, then for periodic perturbations. 

2. Preliminaries 

Let us start with the unperturbed system which is a ring chain T of the form 
sketched in Eig.[Tl choosing the units appropriately we may suppose without 
loss of generality that the circumference of each ring is 27r. We also suppose 
that the particle is exposed to a magnetic field generated by a vector potential 
A; the field is assumed to be perpendicular to the graph plane and homoge- 
neoufl The corresponding vector potential can be thus chosen tangential to 
each ring and constant; since the coordinates we use to parametrize T refer to 

^In fact, the only quantity of importance will be in the following the magnetic flux through 
the rings, hence in general the field must be just invariant with respect to the group of 
discrete shifts along the chain. 
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different orientations in the upper and lower part of the chain, respectively, 
we choose —A as the potential value on the upper halfcircles and A on the 
lower ones. 



Furthermore, in what follows we exclude integer values of A from the 
consideration. This is possible due to the fact that the spectral properties 
of the system are invariant with respect to the change of A by an integer 
which reflects the existence of a simple gauge transformation between such 
cases. We note that in the chosen units the magnetic flux quantum is 27r and 
A = we can then rephrase the above claim saying the systems 

differing by an integer number of flux quanta through each ring are physically 
equivalent. In this sense the case of an integer A is thus equivalent to the 
non-magnetic chain treated in [DET08) . 

The particle Hamiltonian —Aq,,a acts as (—iV — A)^ on each graph 
edge, with the domain consisting of all functions from the Sobolev space 
iFjQc(r) which satisfy the boundary conditions (11.11) at the vertices of T, 
where the quasiderivatives are equal to the sum of the derivative and the 
function value multiplied by the tangential component of eA as usual |KS03) . 
The peculiarity of the present model is that we have two pairs of vectors 
of opposite orientation so their contributions cancel and the left-hand side 
of (jl.ll) is in fact nothing else than the sum of the derivatives taking into 
account different coordinate orientations, in other words, (ffH) is the standard 
S coupling. In contrast to [EM 15] we assume throughout this paper that the 
coupling constant a is the same at each vertex and we are going to determine 
the band-and-gap structure of the spectrum. 



Figure 2. Elementary cell of the periodic system 

In view of the periodicity the band-and-gap structure of the spectrum of 
—^a,A can be computed using Bloch-Floquet decomposition [BK131 Sec. 4.2]. 
Consider an elementary cell with the wave function branches denoted as in 
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Fig- [Hand inspect the spectrum of the Floquet components of — Since 
the operator acts as —:= (i^ — A)^ on the upper halfcircles and as 
— := (i^ + A)^ on the lower ones, each component of the eigenfunction 

with energy F := 7 ^ 0 is a linear combination of the functions 

on the upper graph links and of on the lower ones. In what follows 

we conventionally employ the principal branch of the square root, that is, the 
momentum k is positive for E > 0, while for E negative we put k = in with 
K > 0. For a given E ^ 0, the wave function components on the elementary 
cell are therefore given by 

= e-'^“ {C+e'^^ + C'-e-*'^^) , x G [0, tt) , 

(^(x) =e'^“(i:'+e''=“ + £»-e-‘'^'^) , xG[0,7r). 

As we have said, at the contact point the (5-coupling is assumed, i.e. we have 


V'(-o) = V'(+o) = </5(-o) = </5(-f0), 

— ^'0(—0) -I- ^lp{+0) — ^ip{—0) + ^ip{+0) = Q:'0(O) , 


where the left limits, '(/'(~ 0 ) stc., refer to values on the left-neighbor circle 
parametrized by x G [— tt, 0). On the other hand, the function values at the 
‘endpoints’ of the cell are related the Floquet conditions. 


^(0) = e'®'i/’(^) 1 S>il){—Q) = — 0), 

1^(0) = e'®i^(7r), = e‘®^(/j(7r — 0), 


with 6 running through [—7r,7r); alternatively we may say that the quasimo¬ 
mentum ^6 runs through ^), the Brillouin zone of the problem. In both 
cases the vector potential contributions subtract and S> in (EHl) and (EH) can 
be replaced by the usual derivative. 

Using the explicit structure of the wave function m, its continuity at 
the graph vertices as well as Floquet conditions (EH, one easily deduces that 

C+ = n{A)C- , D+ = fi{-A)D- , (2.4) 


where 


^ gi(0 —Att —feTr) 

^ _ gi(0—A7r+fc7r) 

Thus from the first line in (|2.2p we conclude that 


D- = 


1 -I- fi{A) 

l + K-A) 


(2.5) 


while from the second one we find after simple manipulations the character¬ 
istic equation for —namely 


sinfcTTCOS A 7 r(e^'® — 2^(A:)e‘® -I- l) = 0 (2.6) 


with 

:= -^ ( cos/cTT -I- ^ sinfcTT) . 

^ ^ cos Att V 4A: / 

The special cases A — i G Z and fc G N lead to infinitely degenerate eigen¬ 
values. Assume that none of this possibilities occur, then we get a quadratic 
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equation for the phase factor e'® having real coefficients for any k G ]RUiR\{0} 
and the discriminant 

D = mk?-^)- 

One has to determine the values of for which there is a 0 G [—tt, tt) such 
that (12.61) is satisfied, in other words, those for which it has, as an equation 
for the unknown e‘®, at least one root of modulus one. Note that a pair of 
solutions of (12.61) always gives one when multiplied, regardless the value of k, 
hence either both roots are complex conjugated numbers of modulus one, or 
one is of modulus greater than one and the other has the modulus smaller than 
one. Obviously, the latter situation corresponds to a positive discriminant, 
and the former one to the discriminant less or equal to zero. We shall not 
repeat the discussion which the reader can find in [EM 15) and limit ourselves 
to quoting its conclusions. 

Theorem 2.1. If A — ^ G Z; then the speetrum of —/S.a,A consists of two 
series of infinitely degenerate eigenvalues, namely {k^ G M: f,{k) = 0} and 
{k^ G M: fc G N}. //G K is such that f,{k) = 0, then the corresponding jth 
eigenfunction is given by 

sin kx , ipj+i{x) = smk{7T — x), a;G[0,7r], 

ipj{x) = sinkx , sinA:(a; — tt ) , a;G[0,7r], 

on the jth and {j + l)st circles and vanishes elsewhere for all j G Z. G M 
is such that fc G N, then the corresponding jth eigenfunction is given by 

tpj{x) = sinkx , ifj+i^x) = sinkx , a;G[0,7r], 

'Pj{x) = —sinfci, {pjj^i[x) = (—sinfca;, x G [0,7r], 

on the jth and j + 1st circles and vanishes elsewhere. 

On the other hand, suppose that then the spectrum of —Aa^A 

consists of infinitely degenerate eigenvalues equal to kf with fc G N and with 
the above eigenfunctions, and absolutely continuous spectral bands with the 
following properties: 

Every spectral band except the first one is contained in an interval 
(n^, (n + 1)^) with n G N. The position of the first spectral band depends 
on a, namely, it is included in (0,1) if a > 4(| cos Att) — l)/7r or it is negative 
if a < —4(1 cosAttI + l)/7r, otherwise the first spectral band contains zero. 

We finish the section with calculating (for each circle) the quantum flux 
of particle going through upper and lower semicircles, which are obtained 
from the formulae J.,/, = —i(V^*^^ — ^-nd J(p = 

with wave functions (1). Direct calculations give us 

J^ = 2fc(|C+p-|C-p), 

J^ = 2fc(|D+|2-|D-|2). 

Notice that from the Bloch-Floquet theorem it follows that the above prob¬ 
ability currents are conserved along the chain. On the other hand, exploiting 
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formulae (12.41) . (12.51) . one concludes that 

HA)\-^-l fi{-A) +1 ' 

Jv Im(-^)P-i m(^) + i 

As a result, in the non-magnetic case the probability current has a mirror 
symmetry, J^/J^p = 1. On the other hand, if a nontrivial magnetic field, 
A ^ Z, is applied, the transport becomes in general asymmetric - cf. also 
Rem. I3.2f iil below. 


3. Dual discrete operators 


In this section we consider more general system, namely, we assume that 
the chain graph under consideration has circles of different lengths, and the 
coupling constants as well as the values of the magnetic field could be different 
on different circles. Denote by a = {a^lj-gz and A = arbitrary but 

fixed sequences of real values for coupling constants and magnetic fields; 
both semicircles of the jth graph circle could be associated with the interval 
Ij = {xj,Xj+i) of the length and we require that the following condition 
holds 

inf^>0, (3.1) 

i&i 

where £ = {£_,}_,gz- Consider the corresponding magnetic Laplacian — Aq, a 
on such a chain graph, that is, the operator acting as —= (i^ ± A)^ 
on the jth graph ring with the domain consisting of those functions from 
the appropriate Sobolev space that satisfy the S boundary conditions at the 
graph vertices. 


ijjixj) = ipj{xj) = jjj-i{xj) = (Pj-i{xj ), (3.2) 

Sljilij{xj) + S)jipj{xj) - - ^j-iipj-i{xj) = aji}j{xj ); (3.3) 


the derivatives are naturally meant as one-sided ones. As before the wave- 
function component ipj corresponds to the upper halfcircle of the jth ring 
while ipj stands for the lower one. Our aim in this section is to describe 
a bijective correspondence between the operator —Aq,,.£.a in A 2 (r) and a 
certain operator in £ 2 ( 2 ). We are going to write a difference equation every 
bounded (square summable) solution of which gives rise to a bounded (square 
integrable) solution of the Schrodinger equation corresponding to —Aq, ^ a, 
and vice versa, every bounded (square integrable) solution of the Schrodinger 
equation produces a bounded (square summable) solution of the said differ¬ 
ence equation. This connection plays a crucial role in the following sections, 
where we will study the spectrum of for specifically chosen sets a, 

I and A. 

Denote by (^) the general solution of the Schrodinger equation 


( —Ac(,£,A — fc^) 


\(p{x,k)J 


= 0 , 


"Ak > 0 
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which can be rewritten componentwise on the upper and lower semicircles as 
follows 

/c) = 0 , SJfc > 0 , X €l-i , 

(3.4) 


- k )ip{x, fc) = 0 , 


5JA: > 0 , X . 


Recalling that the wavefunction should satisfy the continuity condition (IQ) , 
we see that the solutions '0 and ip and their quasiderivatives are given by 

ip{x, k) = ^(xj,k) cosk{x — Xj) + 3lj%l){xj+, k) 

^j'ip{x, k) = 0(xj, k) ks\nk[x — Xj) + S>j'ip{xj+, k) cosk{x — xj)) 

(3.5) 

sin/c(a: — xj) 


and 


p{x, k) = (^{xj,k) cos k{x — Xj) + ^jp{xj+, k) 

^jp{x, k) = {—p{xj , k) ksmk{x — Xj) + ^jp{xj+, k) cos k{x — Xj)) 

(3.6) 

for > 0, a: e Xj. Next, let us introduce the vector 


4'j(^,t) = ^ 
and the matrix 

K{k) = 


e'^ip{xj,k) + e '^p{xj,k) 
^j-iipixj — , k) + e~'^S>j-ip{xj — , k) 


' cos k(.j + ^ sin k£j 4 sin kij'' 
/c sin kij + ^ cos kij cos kij ^ 


$,(fc) = 


then taking into account condition dSISl), we conclude that 

K{k)"i!j{k^0) ='i'j^i(k,iAj£j) ^ Sfc > 0, j GZ. 
In a similar spirit we next introduce the vector 

^ ilj{xj,k)+p{xj,k) \ 
^il;{xj-i,k) +p{xj-i,k)J 

to obtain, by virtue of (1531) and (EH), the relations 

cos{Aj-i£j-i)^j{k) = Mj-i{k)'iij{k,0) 

= Lj-i(k)'i>j{k/iAj-i£j-i) , 

where the matrices are defined as follows 


Mj{k) = 
Lj{k) = 


'cos{Aj£j 


0 


, cos k£j — ^ sin k£jj 

^ 1 o'' 

.cos{Aj£j) cos k£j —4 cos(Aj£j) sin k£j 
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Finally, we define the matrix 


Njik) = 


cos{Aj-i£j-i) 


L,{k)K{k){M,_^{k)) 


-1 


cosiAjlj) 

1 /^sin(Hj)-L 


sin(fe^j) cos(Aj_i^j_i)'' 


^3/ ' sin(k£j-i) 

cos{Ajij) 


cos{Aj£j) y 

where fe G ^ := {z: 5>z > 0 A z ^ Z}, to get the relation 

-/Vj(/c)$j(fc) = $j+i(fc), kGA, 


sin(fc£j_i) 

0 


which by continuity of the wavefunction at the graph vertices can be rewritten 
as 


N,{k) 


( i’jik) \ 
Wj-i(fc)/ 


V J ’ 


/c G , 


or equivalently as 


sin(fc£j_i) cos{Aj£j)iljj+i{k) + sin{k£j) cos{Aj-i£j-i)tpj-i{k) 

= (^ sin(fc£j_i) sin(fc£j) + sinfc(£j_i + j k€A, (3.7) 

where tpj{k) := k). Now we can summ up the above calculations. 


Theorem 3.1. Suppose that aj,£j,Aj G K for j G Z and that (EH holds. 
Then any solutions if^x^k) and ip{x,k), k'^ G M., k € A, of (13.41) satisfy 
relation EH- Conversely, any solution of difference equation EH defines 
via 


^(x, k) = tjjj^k) cosk{x — Xj) 


'3 

\Ait 


+ ^ — ifjik) cosk£j) 


sin k{x — Xj) 
sin 


(p{x,k) = ^^fjjik) cosk{x — Xj), 


(3.8) 


+ {'if>j+i{k)e — tpj{k) cos k£j) 


sin k{x — Xj) 
sin k£j 


X Glj, solutions of equations EH) satisfying S-coupling conditions EH, EH- 
In addition, G Lp(T) if and only ifffijfkfj^j^i G ^p(Z) forp G {2,oo}. 


Proof. It remains to prove the last statement (in both directions). Let k^ G K, 
k G A, and assume all the solutions 'ij}{-,k), (p{-,k) and V'j(^) to be real. If 
G Lp(R) and thus I^j^p G Lp(R) we infer ^jif G Lp(R) for all 

1 < p < oo and j G Z. Then {ifjj{k)}j^z G ^p(Z) follows from 

ip{xj,k) = ^^(x, k) cosk{x — Xj) — ^jil){x, k) ^ 
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X €lj, for p = oo, and from 

+ ^{9jip{x+ ,k)f 

= k)f + k)f^ , xGl,, 


for p = 2. Conversely, assume {^j(fc)}jgz G •^p(Z) holds for p = oo or p = 2. 
The case p = oo directly results from (13.81) and the case p = 2 follows from 
(IXa in combination with 


g2iA,(x-x,) ^ ^ k)f 


= i'^iXj,k)f 


AAj-K 


il’j+i — tjjj (k) cos kn 

sin kiT 


X €lj ; 


this completes the proof. 


□ 


Before proceeding further, let us note that using Theorem 13.11 one can 
rewrite formulae for probability current on the jth circle in terms of 
specifically 

2k 

Jipik) = ——— {5i^j(A:)(3fi'!/'j+i(fc) sin AjTT + 3'0j_i_i(fc) cosA^tt) 

Sm fCTT 

— 3^/>j(A:)(5Ji/'j-|-i(A:) cos Ajir — 3i/;j_|_i(fc) sinAljTr)} , 

2k 

Jip{k) = —^sin AjTt — ^tjjj+i{k) cos A^tt) 

sm fCTT 

— cosAjir + sin AjTt)} . 

Remarks 3.2. (i) Since we consider the above formulae for k G A, the corre¬ 

sponding denominators never vanish. On the other hand, integer values 
of k lead to a ‘compact’ sequences {tjjj{k)}j^z with a finite number of 
nonzero elements. 

(ii) In the nonmagnetic case when Aj G Z, the probability currents on upper 
and lower edges with the same vertices are equal, namely, 

Mk) = Mk) = - ^i^j+i{k)<Ai^,{k)) . 

sm /vTT 

Note that this is related to the symmetry of the 5 interaction involved. 
If we replace it by an asymmetric coupling, interesting and [ possibly 
important ‘switching’ patterns between the upper and lower parts of 
the graph may occur [CP 15] . As we have seen in the previous section, 
a nontrivial magnetic field can lead to another asymmetry. 

(iii) In what follows we will demonstrate that for those k which produce £ 2 - 
sequences {^j(fc)}jgz, the latter can be chosen to be real, whence the 
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probability currents read as follows 

Sin kiT 

In other words, since the coefficients 'tl’j{k)4’j+i{f^) decay for a fixed k 
as |j| —> cx), the probability current is ‘circling’ around such localized 
solutions. 


4. Local perturbations of periodic systems 

After these preliminaries, let us pass to our proper topic. Let the Hamiltonian 
of the periodic system —Aa^A be defined as in Sec. [3 and suppose that 
the system suffers compactly supported perturbations of different types. In 
particular, we consider systems with modified values of the magnetic field 
on two adjacent rings, systems with modified values of the coupling constant 
and the magnetic field on a fixed ring, and finally, systems with one ring of 
a different length. Our goal here is to relate spectral properties of perturbed 
and periodic Hamiltonians. 

Let us start with the essential spectrum. We note that if we decouple 
the straight chain by changing the matching condition to separation ones, 
the essential spectrum will not be affected as the resolvents of corresponding 
differ a finite-rank perturbation. This also means that the essential spectrum 
of the halfchains is again the same, just its multiplicity is one instead of two. 
Since we shall consider perturbations of a compact support, we may cut such 
a perturbed chain at two points at both sides of perturbation support. Irre¬ 
spective of the type of the perturbation, the middle part spectrum is discrete, 
hence since cuts are associated with a finite rank, the essential spectrum is 
again the same as for the straight chain. Our aim now is to look into the 
properties of the discrete spectrum of each perturbed operator. 


4.1. Perturbations of the magnetic field 

Let us consider a periodic system with parameters a = {..., a, a,... }, .^ = 
{..., TT, TT,... } and A = {..., A, A,...}, and suppose that on the indicated 
couple of neighboring rings we modify the values of the magnetic field, say, 
to Ai and A 2 . Without loss of generality we may employ a ring numbering 
starting from the chosen pair labeling it by the indices j = I, 2. The perturbed 
Hamiltonian representing the magnetic field A on each ring of the chain 
except the first two and Ai, A 2 on the first and second rings, respectively, 
will be then denoted by — , 4 ^. Our goal in this subsection is to find 

spectral properties of this operator. 
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According to our circle numeration the matrices Nj from the proof of 
Theorem 13.11 take the form 

Ni{k) =( cos ITT I = --—( COS fcTT + — sinfc-TT) , 

^ ^ \ I 0 / cosAitt V 4fc / 

N2{k) = cosAaTT^ ^^{k) = - \ - f COS fcTT + ^ sin fcTT V 

^ ^ V 1 0 / cosA27r V 4fc / 

/2£(fc) -™^\ 

N,{k) j , 


while for j G Z \ {1, 2, 3} one infers that Nj[k) = N{k), where 


N{k) := 




with ^(fc) being the quantity that appeared in (12.61) . As before we have the 
relations 

4>j+i(fc) = , jGZ, 

from which it follows that 


',+4(fc) = (iV(ft))J<i>4(fc), 

j GN, 

(4.1) 

$4(fc) = N3{k)N2{k)Ni{k)^i{k) , 


(4.2) 

',+i(fc) = (iV(fc))^A>i(fc), 

-jGN. 

(4.3) 


It is clear that the asymptotical behavior of the norms of is determined 
by spectral properties of the matrix N. Specifically, let 4>4 be an eigenvector 
of N corresponding to an eigenvalue then \fj.\ < 1 (or \fj,\ > 1, \fj,\ = 1) 
means that ||4>j|| decays exponentially with respect to j > 4 (respectively, it 
is exponentially growing, or independent of j). At the same time if 4>i if an 
eigenvector of N corresponding to an eigenvalue fj, such that |/i| > 1, then 
||<i>j|| decays exponentially with respect to j < 1 (with similar conclusions for 
\fi\ < 1 and \n\ = 1 ). 

By virtue of Theorem 13.11 the wavefunction components on the j-th 
ring are determined by , and thus, in view of dSD and ( 1131 ) . by $4 or 
$1 depending on the sign of j. If $4 has a non-vanishing component related 
to an eigenvalue of N of modulus larger than 1, or $1 has a non-vanishing 
component related to an eigenvalue of modulus less than 1 , then the corre¬ 
sponding coefficients determine neither an eigenfunction nor a generalized 
eigenfunction of —Aai,A 2 - On the other hand, if $4 is an eigenvector, or a 
linear combination of eigenvectors, of the matrix N with modulus less than 
one (respectively, equal to one), and at the same time 4>i is an eigenvector, 
or a linear combination of eigenvectors, of the matrix N with modulus larger 
than one (respectively, equal to one), then the coefficients <I>j determine an 
eigenfunction (respectively, a generalized eigenfunction) and the correspond¬ 
ing energy E belongs to the point (respectively, continuous) spectrum of the 
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operator —/^Ai,A 2 - To perform the spectral analysis of N(k), we employ its 
characteristic polynomial at energy /c^, 

A2-2e(A:)A + l; (4.4) 

it shows that N{k) has an eigenvalue of modulus less than one iff the dis¬ 
criminant of (HU) is positive, i.e. 

ie(fc)l > 1 , 

and a pair of complex conjugated eigenvalues of modulus one iff the above 
quantity is less than or equal to one. In the former case the eigenvalues of 
N{k) are given by 

Ai,2(fc)=e(fc)±v'?(fc)'-1, (4.5) 

satisfying A 2 = A)"^, hence IA 2 I < 1 holds if ^(fc) > 1 and |Ai| < 1 if this 
quantity is < —1. Moreover, the corresponding eigenvectors of N{k) are 

It is convenient to abbreviate 

A*(fc) =C(fc) -sgn(^(fc))V^(fc)2 - I, 

\*[k) = ^(fe) -b sgn(^(fc)) v^^(fc)2 - I. 

The above functions coincide with Ai(fc) and \ 2 {k) if ^(/c) < —1 or with X 2 {k) 
and Ai(/c) if ^(/e) > I, hence the absolute value of A*(fc) is greater than I 
while the absolute value of A*(/c) is less than 1 unless G ct(—A a,^). We 
also introduce the vectors 

=(a 4)) "'“Uw)' 

which play the role of the corresponding eigenvectors. 


Proposition 4.1. Assume that k^ € K. \ (t(—Aq^^),' then k^ is an eigenvalue 
of —Aai,A 2 'Iff the relation 


mxffk) 


[cosAi’nff -\- (cosA27r)^ 
2 (cos Airff 


(4.7) 


is valid for this k. 


Proof. Observe that by virtue of Theorem l3.1l and formulas (14.11) and (14.31) the 
only possibility to construct an eigenfunction of — is by demanding 

that 

$4 ^ It* and $1 ~ M* . 

We conclude from relation (EH) that is an eigenvalue of —Aai,A 2 tff 
det[N 3 {k)N 2 {k)Ni{k)u*{k),uffk)] =0, 
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where the symbol [a, b] here and in the following is used for a 2 x 2 matrix 
with the columns a and b. Taking into account the explicit structure of the 
matrices Nj, we obtain that 


+ -466+ S 

^ cos Bi-n _ cos Att 

cosS27r ^^cosB2'7r 

Exploiting the identities 6 = € co?n ^ = 2^A* — 1 and the fact that 

I A* I > 1, one can easily rewrite the above determinant in the form of (14.7L 
completing thus the proof of the proposition. □ 




Recall that we consider those from the real line that are not in the 
spectrum of —Aa,A- For such values of the energy the function ,^A* is greater 
than one, hence the equation ( 031 ) admits no solution unless 

(cosyli7r)2 + (cosA27r)2 
- 2(cos^7r)2 - 

Let us note that condition (I4.8|l is also sufficient for eigenvalue existence. In 
fact, if it holds, then the right-hand side of ( 031 ) takes values in (1, 
Exploiting the explicit structure of the functions 6 "we see that the interval 
(1, is in the range of ^A* on every spectral gap, since 


i^xnin) 


+ 


(cosA7r)2 IcosAttIY (cosATr)^ 


-1 > 


1 

(cos A-kY 


Thus we get the following claim. 


Theorem 4.2. The spectrum of —Aai,A 2 coincides with the one for —Aa^A 
unless condition (031) holds. On the other hand, if (Ol) holds, then the essen¬ 
tial spectrum of —Aai,A 2 coincides with ct(—A a,^): and moreover, —Aai,A 2 
has precisely one simple eigenvalue in every gap of its essential spectrum. 

As a direct consequence of the theorem is the following 

Corollary 4.3. Suppose that some Aj, say, A 2 , eguals A. The spectrum of 
—A^j := —Aai,a coincides with the one for — A^^^i unless the following 
condition holds 

|cosAi7r| ^ 

IcosAttI 

On the other hand, if this is the case, then the essential spectrum of —Aa^ co¬ 
incides with (t(—A a,A) and, moreover, —A^i^ has precisely one simple eigen¬ 
value in every gap of its essential spectrum. 


4.2. Mixed perturbations 

We again start with a periodic system and suppose that on a certain circle one 
has a new value of the magnetic field as well as the ‘left’ coupling constant, 
denoted by Ai and ai respectively; we introduce ring numbering starting 
from the chosen one. The Hamiltonian representing this mixed perturbation 

















14 


P. Exner and S. Manko 


will be denoted by — Ao,i,Ai and our goal in this subsection is to demonstrate 
spectral properties of — 

As above we employ the recurrent relations = Nj^j with the 

matrices Nj defined as follows 


Ni{k) = 


COS Att '' 
cos Ai TT 


^2Uk) - 
. 1 0 

, 1 0 


N2{k) =1 ) , N,{k) = Nik) 

for j G Z \ {1, 2}. Here 

£.i(k) = - —( cosA:7r+ ^sin/cTr') . 

cosAittV ik J 


Using the same ideas as in the previous subsection we conclude that the 
characteristic equation is of the form det[A 2 (/i:)Ai(fc)M*(fc), u*(/c)] = 0. Thus 
substituting the explicit expression for the matrices Nj and the function 
we find that k^ is an eigenvalue of — Aq^^Ai if and only if 


(4.9) 

holds for this k. Let / stand for the right-hand side of the above relation. 
Then, as k^ varies from the lower end of a gap in —Aq,a to the upper end, 
f{k) is continuous with respect to k and strictly increasing with respect to k^. 
In particular, if | cos Ai7r| > | cos AttI, /(fe) alternately increases from —oo to 
some positive number or from some negative number to oo, starting with the 
increase from — oo in the first gap (the one below the continuum spectrum 
threshold). The sequence of local maxima (that are positive) is increasing with 
respect to the gap number, at the same time, the sequence of local minima 
(that are negative) is decreasing. On the other hand, for | cos Ai7r| < | cos Att |, 
/(fc) alternately increases from —oo to some negative number or from some 
positive number to oo, starting with the increase from —oo in the first gap. 
In this case the sequence of local maxima (that are negative) is decreasing 
with respect to the gap number, while the sequence of local minima (that are 
positive) is increasing. 


Theorem 4.4. For |cosAi7r| > |cosA7r| and ai > a, the operator —Aq^^Ai 
has precisely one simple eigenvalue in every gap of its essential spectrum, ex¬ 
cept possibly a finite number of odd gaps. On the other hand, for \ cos Ai7r| > 

I cos AttI and ai < a, it has precisely one simple eigenvalue in every gap of its 
essential spectrum, except possibly a finite number of even gaps. In particular, 
for sufficiently small a\ — a positive or negative, there is an eigenvalue in 
every gap. For \ cos Ai7r| < | cos AttI and a\ > a, there is precisely one simple 
impurity state in a finite number of even gaps; for \ cos Ai7r| < | cos Att] and 
ai < a, there is precisely one simple impurity state in a finite number of 
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odd gaps. In particular, for sufficiently small ai — a positive or negative, the 
operator —Aai,Ai has no eigenvalues. 


4.3. Perturbations of geometry 


In the final subsection of the section we assume that the periodic system pre¬ 
sented by —Aa^A suffers a geometric perturbation. Specifically, we suppose 
that the first ring of the chain is rescaled in such a way that his length changes 
to 2£, the equal distances between the two vertices being preserved, while the 
other characteristics of the system such as coupling constants and magnetic 
fields remain the same. Denoting for the sake of brevity the corresponding 
Hamiltonian by — A^, we are going to show that its spectral properties de¬ 
pend, in particular, on whether the scaling transformation is a contraction 
{£ < tt) or dilatation {£ > tt), as well as on the sign of the coupling constant. 

Using the same reasoning as in the previous subsections, we first find 
the matrices Nj entering the recurrence relations = Nj^j. For the 

perturbation in question we have 


Ni{k) = 


'26 (fc) - 


sin kl cos Att '' 


i k-K c 


iA£ 


/ 26 (fc) 

N,ik)=i 


sin /cTT cos Al ^ 


1 k£ c 


5 Att 


with 


6(fc) 


1 / sinfc(7r + 6 g ' 

COS \ 2 sin kiT Ak 


Uk) = 


1 / sin fc(7r + £) a 


sin fcTT 


cos Ak \ 2 sin kl Ak 

while for j £ Z \ {1, 2}, we get Nj{k) = N{k). Furthermore, the eigenvalues 
of — A^ are determined by the characteristic equation 

det[A2(fc)Ai(fc)it*(fc), M*(fc)] =0. 

Taking into account the explicit structure of the matrices Nj , we find that 

\2 sin COS Att sin /ctt cos A.^ rt/- / ja sin k£. cos Att \ 


N2ik)N^ik) = 




sin fcTT cos A£ sin k£ cos Att 
sin k£ cos Att 


-26(fc)S 


26(fc)' 


sin k-n cos A£ 
sin k£ cos Att 


' sin/cTT cos A£ sin/cTTCosA^ / 

from which we conclude that € R \ (7(—A„_^) is an eigenvalue of —Ag if 
and only if 

, sin k£ cos Att 


2fi{k) - X,{k)- 


' ■ 1 = 1 - ( 4 - 10 ) 

sm KTT cos A£ 

We restrict ourselves to discussing solutions of equation ()4.10|1 below the 
continuum spectrum threshold, i.e. in the first spectral gap of —Aa^A- As 
varies from the lower end of a gap to the upper end, the left-hand side 
of (|4.10l) is continuous with respect to k. Suppose first that a is negative, then 
the left-hand side of (14.101) is strictly decreasing with respect to and in the 
case of contraction (dilatation) its minimal value is less (respectively, greater) 
than one, hence one obtains one (respectively, no) solution to equation (Id.lOL 
Assume next a > 0. In the case of a contraction the left-hand side of (14.101) 
is strictly decreasing with respect to k^ and its local minimum is greater 
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than one. At the same time, in the case of a dilatation the function first 
strictly decreases to its local minimum the value of which is less than one, 
and then it increase to its local maximum. Moreover, for a fixed natural n 
and a sufficiently large i the function on the left-hand side of (j4.10|) has 
exactly n local minima and maxima with the following properties: all the 
minimum values are zero, while the sequence of maximum values is strictly 
decreasing and greater than one. This means, in particular, that for such an 
^ one obtains 2n eigenvalues in the first spectral gaps. 

Denote by the function counting eigenvalues of — in the first spec¬ 
tral gap of its continuous spectrum. 

Theorem 4.5. For a > 0 and i G (0,7r), we have = 0. On the other hand, 
if a > 0 and i > tt, then > 1, and moreover, ^ oo holds as £ ^ oo. At 
the same time, for a < 0 and £ G (0, tt) we have = 1. If a < 0 and i > tt, 
then = 0. 


5. Weak perturbations of periodic systems 

As before let —stands for the Hamiltonian of the periodic system de¬ 
scribed in Sec. [H Now we are going to discuss situations when the system 
suffers some weak perturbations. Specifically, we suppose that the param¬ 
eters of the system are of the form A + eAj and a + eaj with e G (0,1) 
and ask about the spectrum of the perturbed Hamiltonian in the asymptotic 
regime e —>■ 0. First we focus on perturbations supported on a compact sub- 
domain of the graph and demonstrate the behavior of the eigenvalues in the 
spectral gaps of the periodic operator. Next we turn to systems where the 
perturbation is also periodic and show that in this situation a version of the 
well-known Saxon-Hutner conjecture is valid. 

5.1. Local perturbations of magnetic fields and coupling constants 

Our aim here is to compare spectral properties of — A^ a with those produced 
by a weak finite-rank perturbation. To be specific, we suppose that the cou¬ 
pling constant perturbation strength is eaj, j = 1,..., n, at the vertices with 
the coordinates {tt, ..., nTr}, and at the same time, the magnetic field suffers 
a weak perturbation, namely for the ring indices j from {1,..., n} we have 
an ‘additional’ magnetic potential eAj. The perturbation is controlled by the 
small parameter e and the perturbed Hamiltonian will be denoted by — A^. 
In view of the compact support, the essential spectrum is not affected as one 
can check using the argument used in the opening of Sec. 01 

We are going to demonstrate that, as e —>■ 0-1-, the presence of the 
eigenvalue in the gap of CTess(—A^) is determined by the signs of 
and ^3 S'® of oot Att. With this aim in mind, we mimic the 

argument from the previous section which yields the relation 

4>j+i(fc) = Nj{k,£)<I>j{k), 
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where the matrix Nj takes the form 


Nj{k,e) = 




cos(Aj7r) 

0 


with 




1 


cos(yl + eAj)TT 


cos fcTT + (a + eaj) 


sin fcTT 

Ak 


Note that the above matrix admits the following asymptotic expansion 
Nj = N + ETT tan ATrikf,- + 0{e^) 

as e —>■ 0, where 


N{k) = 


'2m m 


1 


0 


and Mj{k) = 


mm 


ctj sin kn 
27rfc sin Att 


0 


Aj-i — Aj 
0 


for j e Z, we just put Aj = 0 for j G Z \ {1, 2,..., n}. We are going to 
combine these relations with 


^n+2+j{k) = {N{k)y ^ri+2{k) , 

$n+2(fc) = Nn+iik)N„{k). ..Ni{k) $i(fc) , 


j€N, (5.1) 
(5.2) 


N.(k) 

$,+i(fc) = (iV(fc))^$i(fc), -j€N. (5.3) 

It is clear that the asymptotic behavior of the norms of the vectors is 
determined by spectral properties of the matrix N. 

Note that in view of Theorem 13.II and formulae (15.11) and (15.31) the only 
possibility to obtain an eigenfunction of — is by demanding that 


$ 


n+2 


and $1 


thus we conclude from relation (15.21) that is an eigenvalue of — iff 


det 


mk) 


1 


1 


= 0 . 


^X^,{k)J ’ \X*{k)J 

Now we observe that in the limit e —>■ 0 the product A<n behaves as 


(5.4) 


mk) = [N {k)m + eiT tan At: J2iNik)y mi- 3 ik)iNik)m +Om, 

3=0 

and using than condition (15.41) . one gets 

det[A<n{k)u*{k),u^:{k)] = (A*(fc))"“''^ det[u*(fc), it*(fc)] 

n+1 

+ (A*(fc))"’e:7rtan Att ^ det[AIj{k)u*{k), u^,{k)] + 0{e^). 


It is easy to see that det[it*,M*] = A* — A*, and moreover, 
det[M,u*,u*] = X*(^X*Aj + A^A^-i + 


aj sin fcTT 


2kTT sin At: 
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hence that the characteristic determinant takes asymptotically the form 
A* - A. + 2e^(fc)7rtanAl^g + Oie^) 

Finally, observe that the eigenvalues are given via solution to the equation 

eg{k) + 0{e'^) = f{k), (5.5) 


where 


f{k) = 


cot Att 

TT 



9{k) — ^ Aj + 
i=i 


sin kiT 

AkTT^{k) sin Att 


n 

i=i 


Obviously, the sign of function / is determined by the number sgn(— cot Att). 
Furthermore, as varies from the lower end of a gap in a{—Aa^A) to its 
upper end, f{k) is continuous with respect to k and strictly monotonous 
with respect to fe^, and f{k) tends to zero as ^(fe) —>■ 1, the value being 
attained at the spectrum threshold. 

At the same time, the function g has the following properties: g{k) is 
continuous with respect to k as k^ varies from the left infinity to the right 
infinity of the real line, outside a countable set of the second order jumps 
{k G C+ : ^(fc) = 0}. Furthermore, g{k) is strictly increasing (decreasing) with 
respect to if cot Att > 0 (respectively, cot Att aj < 0). Note 

that g{k) holds as ^ —oo, as well as that g{k) = 

any non-negative integer k and that these point are inflection points of the 
function g. Finally, consider some fixed neighborhood of a fc G N, then the 
slope of the function g in this neighborhood tends to zero as N 9 /c oo. 

The listed properties of the functions / and g allow us to state the 
following result describing eigenvalues in the gaps. 


Theorem 5.1. Assume that cot Att > 0- > 0, the operator 

—Ag has no eigenvalues as e ^ 0-1- except in a finite number of even gaps, 
where it can have one eigenvalue per gap. Similarly, for < 0 

—Ag it has no eigenvalues except in a finite number of odd gaps, where it can 
have one eigenvalue per gap. In particular, for a sufficiently large \ ^j\y 

the operator has no eigenvalues, while for a sufficiently small \ 

'^3 ^ 0 (respectively, kas an eigenvalue in the second 

(respectively, the first) gap. 

On the other hand, let cot Att < 0- ^^e opera¬ 

tor —Ag has precisely one simple eigenvalue as e ^ 0 in every gap of its es¬ 
sential spectrum except possibly a finite number of odd gaps. If X)J=i *^3 ^ 0^ 
the operator has precisely one simple eigenvalue in every gap of its essen¬ 
tial spectrum except possibly a finite number of even gaps. In particular, 
for a sufficiently large \ ^^e operator has an eigenvalue in every 

gap, while for a sufficiently small \ YfJj=i ^ (respectively, 

^3 ^0/^ kas no eigenvalue in the first (respectively, the second) gap. 
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The theorem does not cover particular situations when some of the 
perturbations has zero mean. Let us focus on them and suppose first that 
ckj = Oj then the characteristic equation ()5.5I1 reads 



We see that for a sufficiently small positive e the necessary and sufficient 
condition for eigenvalue existence takes the form 



sgn ^ 


(5.6) 


Theorem 5.2. Suppose that ~ 0- spectrum of — coincides 

with the one for — unless condition (1^ holds. If this is the case, then 
the essential spectrum of —A^ coincides with ct(—A a^^), and moreover, — A^ 
has precisely one simple eigenvalue in every gap of its essential spectrum. 

On the other hand, for = 0 the characteristic equation (15.51) reads 


+0(£^) = ~ e-tO. (5.7) 


As kf varies from the lower end of a gap in cr(—A q,_^) to its upper end, the 
right-hand side of this relation is continuous with respect to k and strictly 
increasing with respect to k^. In particular, it alternately increases from —oo 
to zero or from zero to oo, starting with the increase from —oo to zero in the 
first gap, i.e. the one below the continuous spectrum threshold. 

Theorem 5.3. Suppose that = 0- P^r any £ G (0,1) the essential 

spectrum of —Ag coincides with that of —A^^a- Assume that < 0, 

then the operator — A^ has exactly one simple impurity state in every odd gap 
of its essential spectrum for e -A 0-f. If the sum X)j=i positive, then 

it has exactly one simple impurity state in every even gap of its essential 
spectrum as e —?► 0. 

5.2. Weak periodic perturbations 

Let us turn now to perturbations which are periodic along the chain graph. 
Let A and a be sequences with a period p G N, i.e. 



and ask about the effect of the corresponding perturbations on the spectrum. 
Then difference equation (1,3.71) now becomes 
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with 


&\k) = 


1 


cos(A^-^V) 


cos kn - 


(p) 

, sin k-K 
Ak 


k e M., j G Z. 


The relation between the vector whose position indices differing by p can be 
written as 

p 




k G a, 


where ^j(k) has been described in the previous section, while Nj^\k) is 
defined in the same manner as Nj{k) replacing with Since the deter¬ 


minant of Nj is cos(A^_]^7r)/cos(A^-^'^7r), the matrix nj=i 

terminant, and therefore the product of its eigenvalues is one. By the Floquet- 

Bloch theorem $p+i(fc) = e'P'^®$i(A:) and the eigenvalues are given by 

implying 


tr ( (fc) j =2 cosp7r0, 


> 0 , 90 > 0 . 


(5.8) 


Example. The original unperturbed system is, of course, included. Indeed, 
suppose that p = 1, = A and = a, then the characteristic determi¬ 

nant of the corresponding operator reads as 


cos07r = ^(/c), k G A, 


G [-1,1): 


which is what we get from the condition (12.61) unless A — 1 G Z or A: G N. 

In what follows we will assume that the perturbation is weak and put 

= A^^\e) = A + sAj , Aj.^.p = Aj , j G Z , 

(e) = a -I- eaj , aj+p = aj , j G Z, 

denoting by — the corresponding Hamiltonian. Recall from the previous 
section that in this case we have 

m— 1 

Ayp = NP + ETTtanHTT ^ N^Mp.jNP-^-^ + e>(e^) , 
j=o 

which yields, in particular, 

p 

tvAApik) = iT:[N{k)P) + e7rtanH7r^tr(iV(fc)^“^Mj) -|- O(e^) 

1=1 

Given the definition of A (A:) it is straightforward to check that 

/ Up{m) -Gp_i(c(fc))\ 

\Up.,{m) -Up.^m))} ’ 


N[k)P = 
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where Up are Chebyshev’s polynomials of the second kind. Denoting con¬ 
ventionally by Tp Chebyshev’s polynomials of the first kind, and using the 
relation 2Tp{^) = Up{^) — Up-^iC) we arrive at 

trM^(fc) = 2{Tp{^{k)) + £’Kiim{A'K)^{k)Up-i{^{k))g{k)) . 

to state the result, denote by —^e,j the Hamiltonian of the ‘Kronig-Penney- 
type’ chain graph, p = 1, with the parameters a eaj and A + eAj and with 
the corresponding resolvent set p(—Agj). 


Theorem 5.4. Assume that 

p 

Pp = Rn Pi p(-Aej) 

1=1 

is the intersection of all the spectral gaps of the Hamiltonians of the Kronig- 
Penney-type chain graphs indicated above. Then 

PpCpi-Mf^). 

holds for all e small enough. 


Proof. To begin with observe that energies G R in the spectral gaps of all 
‘Kronig-Penney’ Hamiltonians are now simply characterized by 

which is obviously equivalent to 

|^(fe)(l-I-e7rtan(H7r)pj(/c)) I > 1-I-, e —i 0 , jGZ, (5.9) 

where 

a,sinfc 7 r 

+ 4fc^e(fc)sinH^ ■ 

Thus in order to prove the theorem it suffices to show that the validity of 
dSSl) for a fixed k, ’Ak > 0, and all j = 1,... ,p, implies 

I trM^(fc)| > 2 , 


or equivalently 

\Tp{^{k)) + eTrtan{ATT)^{k)Up-i{^{k))g{k)\ > 1-|-C>(e:^), e 0 . 

To this end we will employ explicit expressions of Chebyshev’s polynomials 
of the first and second kind, namely 

and 

[-] 
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From them it follows that the characteristic determinant of the total Hamil¬ 
tonian —can be rewritten as 

Tp{^) + enta.n{An)CUp-i{^)g 

= (-l)*2P"2*"ia(i)6e(i, C) + {{p - 1) mod 2), 

where 

{2i)\{p - i - l)\ 

“W i!(p-l)! ^ 

and 

It is worth mentioning here that the factor a is equal to one for i = 0 and is 
less than or equal to two for i < [|]. At the same time the factor bg can be 
cast into the form 

p—2i 

M.e)= E n ^(1-I-£7rtan(A7r)5j-^)-f 0(£^), £—>•0. 

jl< --<jp-2i n=l 

From what has been already said we finally conclude that 
\Tp{^ik)) + £7rtan(A7r)^(fc)f/p_i(^(fc))5(fc)| 

> 2P-^ l?(fc)(l + £7rtan(A7r)5j(fc))| 

i=i 

p-2 

-2P-^a{l) E Y[\^{k)i^ + £Trtan{ATT)gj^{k))\ 

jl<---<jp-2 n=l 
P 

— ... — ((p — 1) mod 2) E IC(^)(1 + £7rtan(A7r)gj(A:))| -I- 0{£^) 
t=i 

>1-1- 0{e‘^) , £ —>■ 0. 

By that, the proof of the theorem is complete. □ 

One is naturally interested whether the result remains valid also beyond the 
weak-coupling regime in analogy with Theorem III.2.3.6. of [AGHH05] . At 
present, this question remains open. 
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